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The electronic properties of disordered systems at the Anderson metal-insulator transition (MIT) have been 
the subject of intense study for several decades. Thermoelectric properties at the MIT, such as thermopower 
and thermal conductivity, however, have been relatively neglected. Using the recursive Green's function 
method and the Chester-Thellung-Kubo-Greenwood formalism, we calculate numerically the low tempera- 
ture behaviour of all kinetic coefficients Lij . From these we can deduce for example the electrical conduc- 
tivity a and the thermopower S at finite temperatures. Here we present results for the case of completely 
coherent transport in cubic 3D systems. 
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1 Introduction The Anderson model 1 1 1 is widely used to investigate the phenomenon of localisation 
in disordered materials. Especially the possibility of a quantum phase transition driven by disorder from 
an insulating phase, where all states are localised, to a metallic phase with extended states has lead to 
extensive analytical and numerical investigations of the critical properties of this metal-insulator transition 
(MIT)|2|. 

The one-parameter scaling theory plays a crucial role in understanding the MIT |3J- It is based on an 
ansatz interpolating between metallic and insulating regimes |4|. So far, scaling has been demonstrated to 
an astonishing degree by numerical studies of the Anderson model |5 1. However, most studies focused on 
scaling of the localisation length and the conductance at the disorder-driven MIT in the vicinity of the band 
centre. 

In the present paper, we will show numerical results for cubes of volume L 3 which are consistent with 
scaling also at the energy-driven transition across the mobility edge E c near the band tails. In particular, 
we find that at T = the d.c. conductivity a(E) close to E c is well described by the power-law |4][6) 



a(E) 



1 - 



\ E \ < E c ^ 

0, \E\ > E c . 



Moreover, knowing the explicit energy and system size dependence of the conductance g = aL allows 
us to calculate the low-temperature dependence of thermoelectric properties such as thermopower, Lorenz 
number and thermal conductivity in the fully quantum coherent regime. For these, similarly detailed 
theoretical investigations of scaling at the MIT are still lacking. Here we present an approach based on the 
recursive Green's function technique |7|. 
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2 Thermoelectric Transport Properties Thermoelectric transport properties describe the response of 
a material to the presence of a small external electric field and a small temperature gradient VT. E.g., 
the induced (thermo-)electric field E in the Seebeck effect is directly proportional to VT and the constant 
of proportionality is called thermopower S. In general, the dependency of the electric (j c ) and thermal 
currents (/ ) on a electric field and a temperature gradient up to linear order are 1 8 1 



and 



1 



VT 



Ja = — I \e\L2iE-L22 — 



(2) 



Here e is the electron charge. This defines the kinetic coefficients = 1,2). Additionally, for zero 

magnetic field, the Onsager relation states that T 2 i = L12 |8|. Using the definitions of the thermoelectric 
properties 1 8 1 and Eq. 10 one can express a and S in terms of the kinetic coefficients as 



a = L n 



and 



S = 



\e\TLn 



(3) 



Assuming elastic scattering of independent electrons by static impurities or by lattice vibrations one can 
obtain the kinetic coefficients from the Chester-Thellung-Kubo-Greenwood (CTKG) formulation I9l ll0l . 



00 

L l3 {T) = (-!)*+' J a(E)[E - p{T)]^^ 



df(E,p,T) 
dE 



dE 



\ l ,3 



1,21 



(4) 



where p(T) is the chemical potential and f(E, p, T) is the usual Fermi function. 



3 Model and Numerical Method The Anderson model of localisation is based upon a tight-binding 
Hamiltonian in site representation 

w = £ Ei |i)(i|-5>iK>0'h (5) 

where \i) is a localised state at site i and = t are the hopping parameters, restricted as usual to nearest 
neighbours. The on-site potentials £; G [-W/2, W/2] are uniformly distributed random numbers. 

For calculating the energy dependence of the two-point conductance G2 = e 2 /hg2 we use the recursive 
Green's function method 1 7 1 . This method is based on the Kubo-Greenwood formula for the linear response 
regime and non-interacting electrons 1 10]. In order to obtain correct results for purely elastic scattering, 
metallic leads were attached at both ends of the system. In this case the formalism given above is equivalent 
II II to the Landauer-Biittiker formulation HI 21 . To get the conductance g of the disordered region only, we 
have to subtract the contact resistance due to the leads. This gives 1/g = l/g2 — 1 /N. Here N = N(E) is 
the number of propagating channels at the Fermi energy E 1 1 3 1 . An additional complication arises because 
the number of transmitting lead-modes drops to zero outside the lead-energy band, whereas the disordered 
system still has states at these energies. This mismatch gives rise to a lead-induced reduction in g. The 
problem can be overcome by shifting the energy of the disordered region while keeping the Fermi energy 
in the leads in the lead-band centre. This is equivalent to applying a gate voltage to the disordered region 
and sweeping it — a technique similar to MOSFET experiments. 



4 Conductance Scaling and D.C. Conductivity We set the disorder strength to W/t = 12 and impose 
fixed boundary conditions in the transverse direction 1 14 1. For each combination of Fermi energy and 
system size we generate an ensemble of 10000 samples (except for L = 19 and L — 21, where 4000 
and 2000 were generated) and examine the size-dependence of the average and the typical conductance, 
(g) and (In g), respectively. We find that for E-p/t < —8.2 the typical conductance is proportional to the 
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Fig. 1 (top) Conductivity a vs energy computed from (g) /L 
(□), a linear fit with (g) = aL + const. (•) and a scaling 
according to Eq. Q (solid line). The dashed line indicates 
E c /t — — 8.12. Error bars of {g)/L represent the error-of- 
mean obtained from an ensemble average. Also shown is the 
integrand of Eq. J4j for i = j = 1 and three different temper- 
atures (dash-dotted). 

Fig. 2 (left) System size dependence of the 4-point con- 
ductance averages (g) and (In g) for W/t — 12 and Fermi 
energies are given in the legend. Error bars obtained from en- 
semble average. The dashed lines indicate the fit results to 
Eq. and a linear function to (g) and (In g), respectively. 



system size L and the constant of proportionality is negative. This corresponds to an exponential decay 
of the conductance with increasing L and is characteristic for insulating behaviour. For E-p jt being larger 
than —8.05, (<?} is proportional to L. This indicates the metallic regime. 

Following Ref. 1 14 1 we fit the data in the respective regimes to a standard scaling form, i.e. we assume 

{g)=F{ X L 1 ' v ) and (In .9) = F{ X L^), (6) 

where x is a relevant scaling variable. The results for the critical exponent and the mobility edge are given 
in Tabled The obtained values for both types of averages, (g) and (In g), are consistent. The average value 
of v = 1.59±0.18is in accordance with results for conductance scaling at Ep/t = 0.5 and transfer-matrix 
calculations I51 I14I . 

In order to obtain the anticipated power-law form for the conductivity a(E) in the critical regime, we 
need to assume the following scaling law for the conductance, 

(.9) = f{x v L) . (7) 

Then we expand / as a Taylor series up to order riR and \ m terms of the parameter e = (E c — E) /E c up 
to order tor. This procedure gives 

ftfL) = J2 a™ [x v L) m and x (e) = ]T b n e n . (8) 

m— n— 1 

The best fit is determined by minimising the x 2 statistic. Using riR = 3 and tor = 2 we obtain for the 
critical values, v = 1.58 ± 0.18 and E c /t = —8.12 ± 0.03. These values are consistent with our previous 
fits. The to = 1 term aix" in / corresponds to the conductivity close to the MIT. To estimate the quality 
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average E min /t E max /t n R m R 



v Ec/t 



(g) -8.2 -7.4 3 2 

(Ing) -8.8 -7.85 3 2 



1.60 ±0.18 -8.14 ±0.02 
1.58 ±0.06 -8.185 ±0.012 



Table 1 Best fit estimates of the critical exponent and the mobility edge for both averages of g using Eq. The 
system sizes used were in both cases L — 11, 13, 15, 17, 19, 21. 

of this procedure we calculate the conductivity from the slope of a linear fit to (g) throughout the metallic 
regime, and also from the ratio (g) /L. The resulting conductivities are shown in Fig.^ 

5 Thermoelectric Transport Knowing the energy dependence of er it is straightforward to evaluate Eq. 
@ and to compute the temperature dependence of a and S using Eq. yj. To investigate the consequences 
of deviations from the power-law behaviour of a(E) on S(T, E-p) we only use a = a\x v for energies close 
to the MIT {E-p jt < —7.8). Otherwise an interpolation to a obtained from the numerical data is used. The 
results for the thermopower together with the low- and high-temperature expansions using only Eq. Q are 
shown in Fig. [5] Deviations from the pure power-law result 1151 accumulate above k^T/t w 0.5. 

°KT 1 I i I i I i rj , e c -e f =-o.io 



In Fig.^one can see that at this temperature the integrand in Eq. (0} has already a significant contribution 
from the pure numerical data. Nevertheless, the deviations of S from the expansions are small even for 
higher temperatures. This indicates a certain robustness of the results in Ref. 1151 . 
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Fig. 3 Thermopower S vs a dimensionless temperature 
parameter k^,T /t. Empty symbols indicate metallic, filled 
symbols critical and stars insulating behaviour. The dis- 
tance from E c is given in the legend. Full lines are 
obtained from a high-temperature expansion 1151 . The 
dashed lines indicate S using the Sommerfeld expansion 
1151 . In both cases v and E c obtained from the fit to Eq.0 
are used. 
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